Toward a description of contact line motion at higher capillary numbers 
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The surface of a liquid near a moving contact line is highly curved owing to diverging viscous forces. 
Thus, microscopic physics must be invoked at the contact line and matched to the hydrodynamic 
solution farther away. This matching has already been done for a variety of models, but always 
assuming the limit of vanishing speed. This excludes phenomena of the greatest current interest, 
in particular the stability of contact lines. Here we extend perturbation theory to arbitrary order 
and compute finite speed corrections to existing results. We also investigate the impact of the 
microscopic physics on the large-scale shape of the interface. 



The moving contact line problem is a famous example 
of hydrodynamics failing to describe a macroscopic flow 
phenomenon. But it was only in 1971 that Huh and 
Scriven discovered that the viscous dissipation in the 
fluid wedge bordered by a solid and a fluid-gas interface is 
logarithmically infinite if the standard no-slip boundary 
condition is applied at the solid surface. Thus infinite 
force would be required to submerge a solid body, and a 
drop could never spread on a table. 

This result is of course contradicted by observation, 
and physical effects that relieve the singularity have to be 
invoked near the contact line, which go beyond the stan- 
dard description. A great variety of possible mechanisms 
have been proposed, and indeed there is no reason to be- 
lieve that for different solid-fluid-gas systems always the 
same mechanism is involved. However, a question rarely 
considered is whether the choice of different microscopic 
mechanisms would make a great difference when looked 
at macroscopically. 

In a recent paper Q we compared various microscopic 
models in the case of perfect wetting. We found that the 
length scale that appears in the expression for the inter- 
face shape is strongly speed dependent, in a fashion that 
depends on the model. Here we are going to show that 
this dependence is much weaker in the case of partial 
wetting, and differences only come in at higher order in 
an expansion in capillary number Ca = Ur]/j, where i] 
is the viscosity of the fluid, and 7 the surface tension be- 
tween liquid and gas. Finite capillary number corrections 
are of interest for various situations of "forced" wetting, 
in which Ca is no longer asymptotically small, and previ- 
ous theories for the dynamic interface angle break down 

a 

For simplicity, we perform our calculations within the 
framework of lubrication theory, assuming a parabolic 
flow profile. This limits applications to the case of small 
contact angles, but without altering the essential struc- 
ture of the problem. We consider the neighborhood of 
a contact line moving with speed U across a solid in a 
frame of reference in which the contact line is fixed at 
the origin of the coordinate system (see FigQ. To re- 
lieve the corner singularity, we allow the fluid to slide 



FIG. 1: A schematic of the interface near the contact line: In 
the frame of reference in which the contact line is stationary, 
the solid moves to the right in the "wetting" situation con- 
sidered here. At the contact line, h(0) = 0, the slope of the 
interface is 8 e . 



across the solid surface, following the generalized Navier 
slip law 



u|y=0 — U = A h -7^— |y=0 



dy 



(1) 



at the plate, where h(x) is the thickness of the fluid layer, 
and A is taken as a constant rather than a speed depen- 
dent quantity. The case a — 1 corresponds to the usual 
Navier slip, and Q is a simple generalization involving 
only a single length scale A. The resulting lubrication 
equation is |(j 



3Ca 



h 2 + 3X 2 - a h a 



(2) 



The left-hand side corresponds to viscous forces, di- 
verging as the contact line position h(0) = is ap- 
proached, but weakened by the presence of slip. Vis- 
cous forces are balanced by surface tension forces on the 
right, resulting in a highly curved interface near the con- 
tact line. In comparison, other forces like gravity have 
been neglected. This restricts the validity of @ to a 
distance from the contact line below the capillary length 

= ylKpg)- We also assume that the angle at the 
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contact line h'(0) = 9 e is constant, independent of speed. 
Hence it has to coincide with the equilibrium contact an- 
gle, in order to give the right result at vanishing speed. 

Since we want to investigate the neighborhood of the 
contact line, it is convenient to introduce the scaled vari- 
ables 

h(x) = 3V( 2 -«)AF(0, £ = arf e /[3 1/(2 - a) A], (3) 
which leads to 



H 2 + H a 



-H' 



(4) 



From the scaling Q it is evident that the curvature 
of the interface h"{x) scales like A" 1 , where A is in the 
order of nanometers Thus, in order to match the 
local solution near the contact line to an outer profile 
with a curvature of order l/£ c , the curvature H"(£) has 
to vanish for large £. This means the boundary conditions 
for the solution of J3J are 



H(0) = 0, H'(0) = 1, H"(oo) = 0. 



(5) 



The only parameter appearing in the problem is now the 
rescaled capillary number 6 = SCa/d^. 

For 5 > equations and (0 have a unique asymp- 
totic solution, due to Voinov for which the slope be- 
haves like H'(g) = [351n(£/£ )] 1/3 for £ > 1. This solu- 
tion has vanishing curvature at infinity and only contains 
a single free parameter £ 0l to be determined by matching 
to the contact line. In the present paper, we are going 
to deal exclusively with this wetting situation. If <5 < 0, 
the mathematical structure of Q changes completely. 
This can be seen from considering the simpler equation 
S/H 2 = —H'", valid for large H. Namely, it follows from 
an exact solution 0] to this equation, that for S < 
all solutions have strictly positive curvature at infinity. 
The consequences of this observation for the stability of 
contact lines are explored in another paper [ic| . 

The mathematical problem to be tackled in this letter 
consists in computing; £n as a function of 5. As was done 
m previous works dEllil El, we proceed by expanding 
around the trivial solution at zero speed 6 — 0. In this 
case, the solution of J3J and 101 is evidently given by 
= £. Hence the perturbation expansion we seek 
looks like 



(0) 



for the slope. 

This is to be c omp ared to the full asymptotic expansion 
of m m (£) 14], the leading term corresponding to 
Voinov's solution: 



ff'(0 = [3-51n(£/£o)] 1/3 1 + E 



All coefficients b\ = 1/3, 62 = —10/27,... are readily 
computable. Expanding Q in 8 and comparing to © 
leads to the following structure of the expansion of ln(£ ): 



-j 00 
31n(£o) = T+5>+^ 



(8) 



i=0 



Substituting this back into (J7J), the large-£ behavior of 
the Hi(£) in JBJ is given in terms of the coefficients c$: 

= MO + 1 + C1A 1 

ff 2 (0 = -ln 2 (£)-(2 + 2 Cl /3)ln(0+ k -> 00 (9) 
+C2/3 - c 2 /9 - 2ci/3 - 10/3. J 

To compute c\ we have to solve (0J to first order in 5: 

1 



H'{' 



e 2 + e 



(10) 



with boundary conditions -ffi(O) = 0,H[(0) — 1, and 
H"(oo) = 0. According to |[?JJ, the constant c\ is given 
by 



l + c x /3 = lim ff((£)-ln(£+l). 

£ — >QC 



(11) 



Integrating (|10|) twice, we can thus write 



l + ci/3 = 



OO /•£ 



Joo 
00 /»£ " 



.70 

00 



-(0 



1 



(c+i) 2 
1 

(e+Tj 1 



e- a + i (e + 1) 2 
1 



2 - a 



ln(£ 2 - tt + l)-ln(l + £)- 



i + e 



d£d£ = 
= 1 



independent of a. Remarkably ci, the first-order correc- 
tion to £oi always vanishes, regardless of how the length 
scale A is introduced near the contact line. 

The problem at second order can be tackled in precisely 
the same manner, using the equation for the second-order 
problem 



fl*" = H x (0 



2g + l 
C 2 (C + 1) 



(12) 



where we have specialized to the standard case a = 1 for 
simplicity. The case a — can be treated in exactly the 
same manner. Integrating the first order equation 110|) 
thrice gives 

Hi = [ln(£ + + l) 2 - £ - e ln(0]/2, 

thus specifying the r.h.s. of equation (|12(l . The trick 
used to calculate c\ at first order can be repeated at the 
next order, using the second equation of © and c% = 0. 
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FIG. 2: A comparison of 1(14(1 with the numerical result for 
the characteristic length L, using numerical integration of 
for a — 1. 



Simplifying the resulting double integral as before, we 
find 



c 2 /3 - 10/3 = lim + ln 2 (£ + 1 

£ — 

-6ffi(Q(2e + l) + 2e 3 ln(g + l) 

£ 2 (£ + i) 2 



) + 21n(£ + l) = 
d£ = tt 2 /6 - 7/2. 



In summary, we thus have c-2 = (tt 2 — l)/2 for a = 1 
and C2 = (37r 2 — 4)/8 for a = 0, which follows from 
a very similar calculation. It is evident that the same 
procedure can be repeated at arbitrary order, although 
the calculation rapidly becomes analytically intractable. 

Rewriting the rescaled solution H(£) in terms of the 
physical profile h(x), we find the slope of the interface to 
be 



h' 3 (x) -Q\ = 9Ca\n(x/L) 



(13) 



for any capillary number. This is the form originally pro- 
posed by Voinov [8|, using more qualitative arguments. 
The length L appearing inside the logarithm can be writ- 
ten as an expansion in the capillary number: 



L 



3 l/(2-a) A 



(14) 



Integrating equation @ numerically, a comparison 
with 1(131) can be made, giving L as function of capil- 
lary number. In Fig[3 it is clearly seen that our expan- 
sion 114(1 describes the initial departure from the leading 
order result quite well. However, when the corrections 
amount to about 10% of the leading order, higher order 



terms become important. Thus as a rough estimate, the 
present approach can be trusted if <5 = Ca/Sl < 0.05. 
It would be interesting to systematically investigate the 
dependence of L on the capillary number beyond that 
value. 



Comparing (|T3")) to de Gennes' result ^ 



{h' 2 (x) - 6l)h'{x) = 6Caln(x/A), 



(15) 



one finds that the two laws agree if the departure of h'(x) 
from 9 e is small. Beyond the leading order expansion of 
h'(x) in Ca, however, 1(13(1 and ((15(1 are inconsistent. This 
casts doubts on the original argument for 1(15(1 . which 
should also apply to the class of simple slip models con- 
sidered here, and which has already been reviewed in a 
critical light in . In particular, this calls into question 
de Gennes' theory 0] of contact line instability, which 
crucially uses Q15|). 

Finally, it should be realized that we have assumed 
that the local profile as described by © can be treated 
independently. The only effect of the outer profile comes 
in with the boundary condition of vanishing curvature 
(0. Earlier work [(J [llj includes the matching to an 
outer profile, but arrives at the same results as we do 
to leading order in Ca. At elevated capillary numbers, 
however, different parts of the solution may interact in a 
non-trivial way, and matching may be necessary. 

Thanks are due to Howard Stone for many discussions 
on the subject of contact lines, and to Rich Kerswell for 
help with the analytical calculation of integrals. 
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